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Abstract 

The mapping transformation technique is applied to obtain exact results for 
the spin-1/2 and spin-,!? (S = 1/2, 1) Ising-Heisenberg antiferromagnetic chain in a 
presence of external magnetic held. Within this scheme, a held-induced first-order 
metamagnetic phase transition resulting in multiplateau magnetization curves, is 
investigated in detail. It is found that the scenario of the plateau formation depends 
fundamentally on the ratio between Ising and Heisenberg interaction parameters, 
as well as on the XXZ Heisenberg exchange anisotropy strength. 
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1 Introduction 

Quantum antiferromagnetism in lower dimensions is one of the most fascinating subjects 
in the condensed matter physics. In particular, the antiferromagnetic quantum Heisenberg 
chains (AFQHC) with small spins have attracted much attention on account of the rich 
quantum behaviour they display. Nevertheless, due to the strong quantum fluctuations the 
classical Neel state is not more an eigenstate of the Hamiltonian and thus, more interesting 
quantum phases should be expected to occur in the ground state. The nature of these 
phases, however, basically depends on the spin value of atoms. In fact, as conjectured 
Haldane [1] in 1983, the integer spin AFQHC have a finite energy gap between the ground 
state and the first excited state, while the half-odd integer ones possess a gapless excitation 
spectrum. 

Another striking feature of the AFQHC is the appearance of fractional magnetization 
plateaus in the magnetization process. Extending the original Lieb-Schultz-Mattis theo- 
rem [2] Oshikawa, Yamanaka and Affleck (OYA) [3] argued that the magnetization per 
site m can be topologically quantized as: p(S u — m) = integer, where p is a period of 
ground state in the thermodynamic limit and S u denotes a total spin of an elementary 
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unit. However, this condition represents just the necessary condition for the plateau-state 
formation and does not directly prove its existence. It is therefore of interest to investi- 
gate how the plateau state is related to the periodicity of specific model [4]. Moreover, 
from the theoretical viewpoint the plateau state can also be regarded as a spin-gap state. 
Hence, the zero-temperature magnetization curves with plateaus bring an insight into the 
ground-state properties of the system, since the field-induced spin gaps reflect the gapped 
excitation spectrum. 

Despite of an extensive theoretical effort focused on AFQHC, there still exist only 
few exactly solvable models with pure Heisenberg exchange interactions [5], especially for 
mixed-spin chains [6] . On the other hand, an exact solution for the chains with alternat- 
ing Ising- and Heisenberg-type exchange interactions can be attained in less sophisticated 
manner. Indeed, the exact solution for the Ising-Heisenberg bond alternating chain (orig- 
inally proposed and solved by Lieb et al [2]), has been recently successfully generalized 
to the case of anisotropic Heisenberg interaction [7]. In order to avoid mathematical 
complexities connected with the noncommutability of relevant spin operators, we will in- 
troduce in this article another class of the Ising-Heisenberg chains (with period p — 3), 
that can be treated exactly within the mapping transformation method. However, the 
considered model naturally enables to analyse in detail the mutual competition between 
Ising- and Heisenberg-type interactions and moreover, it also proves to be very useful in 
view of the confirmation of multiplateau magnetization curves by an exact calculation. 

This paper is organized as follows. In Sec. 2 the detailed description of the model, as 
well as the fundamental aspects of transformation technique are presented. In Sec. 3 we 
are concerned with the analysis of the most interesting numerical results for typical spin 
cases and finally, in Sec. 4 some concluding remarks are also drawn. 

2 Model and method 

In this article we will study a mixed spin-1/2 and spin-S* (S = 1/2, 1) Ising-Heisenberg 
chain in a presence of external magnetic field. The structure of the considered mixed-spin 
chain is depicted in Fig.l, where the black circles denote the spin-1/2 atoms and the grey 
ones represent the spin-S* atoms. The total Hamiltonian of the system is given by: 

H = J £[A(Sf + SfS]) + S?S?\ + Ji E Slti -Hj^K-HbY, s z k , (1) 

i,j k,l I k 

where fif and S% (a = x,y,z) denote the well-known components of standard spin-1/2 
and spin-S* (S = 1/2, 1) operators, respectively. The parameter J stands for the Heisen- 
berg interaction between nearest-neighbouring spin-S* atoms (the grey atoms) and A 
is the anisotropy parameter that allows to control the anisotropic XXZ interaction be- 
tween an easy-axis regime (A < 1) and an easy-plane regime (A > 1). Furthermore, 
the interaction parameter J\ describes the Ising-type exchange interaction between pairs 
of nearest-neighbouring spin-1/2 and spin-S* atoms and finally, the terms incorporating 
Ha and H B respectively, describe the coupling of spin-1/2 atoms and spin-S* atoms to 
an external magnetic field. As we can see from Fig.l, all pairs of nearest-neighbouring 
Heisenberg atoms are surrounded by Ising-type atoms only, thus the model under inves- 
tigation can also be viewed as an Ising model the bonds of which are doubly decorated 
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by the Heisenberg atoms. In view of further manipulations, it is useful to rewrite the 
total Hamiltonian 7i as a sum of the bond Hamiltonians, i.e. Tt = Y,k=i^-k, where N 
denotes the total number of Ising-type atoms and the summation runs over all bonds of 
the original (undecorated) chain. The bond Hamiltonian Tik contains all the interaction 
terms associated with the fcth couple of Heisenberg atoms (see Fig.l), and it is given by 

T~tk = J[A(S kl S k2 + -Sfci'S'L) + ^1^2] + ■M'SfciAfci + ^jwAjw) 

-H B (S z kl + S z k2 ) - H A (Ji' kl + AL)/2- (2) 

Now, exploiting the usual commutation relation for the bond Hamiltonians (i.e. [Hk, 7~£j] = 
0, for k 7^ j), the partition function of the system can be partially factorized, namely, 

N 

2 = Tr M II ^s k Jr Sk2 exp(-FH k ). (3) 
fc=i 

In above, (3 = (fcsT) -1 , k B being Boltzmann constant and T the absolute temperature. 
Tr^} means a trace over the degrees of freedom of the Ising spins and Tr5 fel Tr 5fe2 denotes 
a trace over the kth. couple of Heisenberg spins. At this stage one easily observes that 
the structure of relation (3) implies a possibility to introduce the decoration-iteration 
mapping transformation [8] 

Tr 5fcl Tr 5fc2 exp(-/374) = A exp^Rfi^ + (3H (^ kl + ^ 2 )/2]. (4) 

As usual, the unknown transformation parameters A, R and H can be attained by taking 
into account remaining degrees of freedom of both Ising spins (/iki and ^2). In this way 
one obtains for the transformation parameters A, R and H the following expressions 

A=(V 1 V 2 VZ) 1 '\ /3i? = m(^), (3H = PH A -In(^-), (5) 

where the functions Vi, Vi and V3 depend on the spin value of Heisenberg atoms, as well 
as on the parameters of Hamiltonian (1) and they are summarized for both investigated 
spin cases in the Appendix. 

Here, one should emphasize that the mapping relations (4)-(5) enable to transform 
the Ising- Heisenberg mixed-spin chain onto a simple spin-1/2 Ising chain with an effective 
exchange parameter R, placed in an external magnetic field of the magnitude H . Indeed, 
substituting (4) into (3) gives the following equality 

Z((3, J, J ls A, H A , H B ) = A N Z (P, R, H ), (6) 

which relates the partition function of Ising-Heisenberg chain Z and that one of the 
spin-1/2 Ising chain Z . Since the explicit expression for Z is well known [9], we can 
then straightforwardly calculate all relevant thermodynamic quantities. For example, the 
Gibbs free energy Q of the mixed-spin Ising-Heisenberg chain is given by 

Q — Qq — Nk B T In A, (7) 

where Q§ = — k B T In Z denotes the Gibbs free energy of the corresponding spin-1/2 
Ising chain. Next, by differentiating the Gibbs free energy Q with respect to H A and 
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Hb respectively, one directly obtains solution for the total sublattice magnetization. Of 
course, other thermodynamic quantities can also be calculated on the basis of familiar 
thermodynamic relations, e. g. the entropy S and the specific heat C can be calculated 
from 

*=-(§f)-. c = - r &- < 8 > 

Nevertheless, a similar thermodynamic approach cannot be used for the calculation of 
other important quantities such as staggered magnetization, quadrupolar momentum or 
some correlations. Fortunately, Eq. (6) in conjunction with the transformation formula 
(4) allows after an elementary algebra the derivation of following exact spin identities [10] 

1 f ( oa o7 ;2 ;2 \\ _ /^Sk^SwhiSkl-i ^k2i Afcl) Ajw) ex P( — PUk) \ / n \ 

b k2 , fan ^2)) - \ — ; — /> l y J 

x TTs kl Tr Sk2 exp{-(3H k ) ' 

with arbitrary function f\ depending exclusively on Ising spin variables and the function f 2 
depending on the spin variables from the kth bond only. The superscript a, 7 = (x, y, z) 
label the spatial components of spin operators and finally, the symbols (...) and (...) 
stand for the standard ensemble average in the Ising-Heisenberg and its equivalent simple 
Ising model, respectively. However, the above spin identities considerably simplify the 
calculation of a large number of quantities. Indeed, for the reduced sublattice magne- 
tization (m z A , m z B ), the total single-site magnetization m and the staggered sublattice 
magnetization (m A , m s B ), one attains after straightforward algebra 

1 1 

m A = 2^*1 + = 2^*i + ^2)0 = m , 

m% = \{S z kl + S z k2 ) = (V./V, - V 5 /V 2 + 2V 6 /V 3 )/2 

-2m (^M + V 5 /V 2 ) + 2e (V A /V 1 - V 5 /V 2 - 2V 6 /V 3 ), 

m = (m A + 2m z B )/3, 
1 1 
m A = gxAfci - ^2) = ~ ^2)0 = mo. 

\(S z k , - SQ = -m s V 7 /V 3 . (10) 



In above, mo, an d £ o represent the reduced magnetization, staggered magnetization 
and nearest-neighbour correlation of the corresponding undecorated Ising chain and the 
coefficients V1-V7 are listed for both investigated spin cases in the Appendix. 

Finally, let us define some pair correlation functions and the quadrupolar momentum, 
which are also very useful for understanding of magnetic properties of the system, namely, 

Ihh = (Skl^te) = ('S'fcl'S'L)' Qhh = (^1^2)) Qu = (AfclAjw) = e o> 

1th = + S z k2 fi z k2 ), V = \{{S z kl f + {S z k2 f). (11) 

In these equations, the subscripts denote the type of atoms and superscripts the space 
direction. One should also notice that the definition of the parameter 77 is obviously 
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meaningful for S > 1 only. Although, the derivation of relevant equations for these 
quantities is straightforward, the calculation procedure by itself is rather lengthy and 
tedious, thus the details are not presented here. 

3 Numerical results and discussion 

Before discussing the most interesting numerical results, it is worth mentioning that some 
preliminary results for the ferromagnetic version of the model (J < 0, J\ < 0) have 
already been published by the present authors elsewhere [11]. For this reason, in this 
article we will restrict our attention to the doubly decorated Ising-Heisenberg chain with 
antiferromagnetic interactions only (i.e. J > 0, J\ > 0). The particular attention is 
focused on the ground-state analysis and the appearance of plateaus in the chains with 
different decorating spin S. Among other matters, we will directly prove the existence 
of double-plateau magnetization curve in the spin S — 1/2 chain, more precisely, the 
magnetization curve with plateaus at m = and 1/6. On the other hand, in the spin 
5=1 chain a greater diversity of magnetization process will be confirmed, in fact, we 
will prove the existence of double-plateau (m — and 1/2), triple-plateau (m = 0, 1/6 
and 1/2) and quadruple-plateau (m = 0, 1/6, 1/3 and 1/2) magnetization curves. 

In addition, since each couple of the Heisenberg atoms is surrounded by the Ising 
atoms only, the relevant spin deviations cannot propagate through the Ising bonds and 
therefore, the quantum fluctuations are necessarily localized within the unit cell (within 
the four-spin cluster consisting of the Heisenberg spin pair and its nearest-neighbouring 
Ising spins). Owing to this fact, the observed plateaus cannot be considered as the OYA 
plateaus, i. e. the plateaus with collective eigenstate extending over the whole chain. 
Nevertheless, it is interesting to note that all observed fractional magnetization satisfy 
the OYA condition for p = 6 (the period of translational symmetry should be twice as 
large as the periodicity of Hamiltonian as a consequence of an antiferromagnetic nature 
of ground state). 

3.1 Spin S = 1/2 chain 

We begin our analysis by considering the effect of exchange anisotropy A and uniform 
magnetic field (i.e. Ha = Hb = H) on the ground-state phase boundaries of the spin 
S — 1/2 chain. For this purpose, we have displayed in Fig.2 some typical ground-state 
phase diagrams in the A — H/J plane for J\j J = 1.0 and 2.0. As one can see from 
this figure, the relevant phase boundaries separate three or four distinct phases, namely, 
the antiferromagnetic (AF), ferrimagnetic I (FI), ferrimagnetic II (FII) and saturated 
paramagnetic (SP) phase. One also observes that both ferrimagnetic phases (FI and FII) 
represent an intermediate phase between AF and SP phase occurring due to the first-order 
metamagnetic transition. As we have already mentioned before, different phases can be 
distinguished by analysing the magnetization and the correlation functions at T = 0. In 
this way one finds the following ground-state results for particular phases: 
Antiferromagnetic phase (AF): 

(q x h lQ Z h h ,Q Z ih,Qu,m s A ,m s B ) = (-JAQ~\ -1/4, -hQ~\ -1/4, 1/2, -2J.Q- 1 ); 
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where we have defined the function Q = 4y J 2 + (JA) 2 , in order to write the relevant 
expressions in more abbreviated and elegant form. 
Ferrimagnetic phase I (FI): 

(q x h x h ,q z h z h ,qH,qti Z ,m z A ,m z B ,m) = (-1/4,-1/4,0,1/4,1/2,0,1/6); 
Ferrimagnetic phase II (FII): 

(<ffi,&,<fi,&,rn' A ,m%,m) = (0,1/4,-1/4,1/4,-1/2,1/2,1/6); 
Saturated paramagnetic phase (SP): 

teC^,C^,m^,m) = (0,1/4,1/4,1/4,1/2,1/2,1/2). 

Moreover, it is also noteworthy that all ground-state phase boundaries can be expressed 
analytically as follows (see Fig. 2), 

a) for Jx/J < 1.0 

(1) H/J = tJA 2 + {JJjy - A, (2) H/J = (A + J x / J)/2 + 1/2, (12) 

b) JJJ > 1.0 

(1) H/J = ^A 2 + (Jx/J) 2 - A, (2) H/J = (A + J l /J)/2 + 1/2, (13) 

(3) H/J = ^A 2 + ( Ji/ J) 2 - Ji/J + 1, (4) H/J = JJ J, (5)A = Jx/J-l. 

In order to demonstrate the diversity of the magnetization process, we have plotted in 
Fig. 3 the low-temperature magnetization curves for various exchange anisotropies A. The 
detailed examination of these dependencies reveals that the mechanism of the magneti- 
zation process depends basically on the ratio between Ising and Heisenberg interaction 
constants, more specifically, on whether J 1 /J < 1.0 or J 1 /J > 1.0. In the former case, the 
plateau state arises due to the alignment of the Ising spins towards the external-field di- 
rection (FI phase) regardless of the anisotropy strength A (see Fig. 3a). Naturally, in this 
case there is no other possibility for the formation of intermediate plateau. On the other 
hand, in the case of J\/ J > 1.0 the metamagnetic transition to the FI phase can be ob- 
served for the stronger anisotropies A only (see Fig.3d), while for the weaker anisotropies 
the FII phase becomes the stable one. In this phase, the Heisenberg (Ising) spins align 
parallel (antiparallel) with respect to the external-field direction as it is apparent from 
Fig. 3b. Moreover, the situation for the most interesting point at which both intermediate 
phases (FI and FII) coexist, is depicted in Fig. 3c. Referring to this plot, the coexistence 
of both ferrimagnetic phases is also reflected in the mixed feature of the magnetization 
curve (compare magnetization curve from Fig.3c with those in Figs. 3b and 3d). 

In order to enable an independent check of the magnetization scenario, we have also 
studied the relevant low-temperature dependencies of the nearest-neighbour correlation 
functions introduced in Eq. (11). In Fig.4a, the variations of the correlation functions 
with anisotropy A are shown for the system without external magnetic field. As one 
can see, the correlation function qff h between Heisenberg spins takes its saturation value 
irrespectively of A, what means that all nearest-neighbouring Heisenberg spin pairs align 
antiparallel with respect to each other. Moreover, the saturated value of correlation 
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qff = —1/4 indicates a perfect antiferro magnetic alignment also in the Ising sublattice 
(between third nearest-neighbour Ising spins). Contrary to this behaviour, the perfect 
antiparallel alignment between Ising and Heisenberg spins is destroyed as A increases 
from zero (see the correlation q^). Hence, the Ising and Heisenberg spins are oriented 
randomly with respect to each other, the degree of randomness being the greater, the 
stronger the exchange anisotropy A. In addition, it is clear that the anisotropy term A 
is also responsible for the onset of an interesting short-range ordering (nonzero qf£) in xy 
plane. Anyway, the value of correlation qf^ can be thought as a measure of the strength 
of local quantum fluctuations appearing in the spin system. Furthermore, in Figs. 4b- 
4d we have displayed the field dependencies of the correlation functions corresponding 
to the magnetization curves from Figs.3b-3d. The depicted behaviour for the correlation 
functions is in complete agreement with results for the magnetization curves and moreover, 
it enables better understanding of the magnetic ordering in relevant phases. As a typical 
example we can mention the ferrimagnetic phases. Although the total magnetization 
of both the FI and FII phases is the same, the results for other quantities indicate a 
fundamental difference between them. In fact, in the FI phase the pairs of Heisenberg 
spins create singlet dimers and thus, they do not contribute to the total magnetization 
which is nonzero due to the fully polarized Ising spins only, as already stated before. This 
observation would suggest that the formation of FI plateau state should be based on the 
quantum mechanism with a significant influence of local quantum fluctuations. In contrast 
to this, the FII phase represents the standard ferrimagnetically ordered phase usually 
observed in the pure Ising systems. Indeed, the generation of the FII plateau state is 
nothing but the gapped excitation from the Neel state, implying the magnetization process 
with a 'classical' Ising-like mechanism. Finally, one should notice that after exceeding the 
saturation field given by conditions (12) and (13), the ground state becomes fully polarized 
(SP phase) and all spins are aligned in the external-field direction. 

Now, let us investigate the finite-temperature behaviour of the system. Firstly, we 
take a closer look at the thermal dependencies of total magnetization that are shown in 
Fig. 5. As one can see, the initial value of total magnetization takes one of three possible 
values m = 0, 1/6 or 1/2 for the AF, FI (FII) or SP phase, respectively. Furthermore, 
there are two special cases which correspond to the coexistence of the relevant phases. 
In these cases, we have obviously obtained m = l/12orl/3atT = 0. It is also easy to 
observe here that the most interesting dependencies appear for external fields from the 
neighbourhood of phase boundaries. Hence, the observed rapid increase (decrease) of the 
magnetization can be attributed to the thermal excitations of huge number of spins which 
occur due to the competitive influence of both phases. 

For completeness, we have also plotted the entropy (Fig. 6a) and the specific heat 
(Fig.6b-6d) as a function of temperature. As one can expect, the entropy of the system 
does not vanish for the boundary external-field values H/J = 1.0 and 2.25, at which the 
relevant phases coexist in the ground state (see Fig.6a). However, we should mention 
that for any other external fields the entropy vanishes as the temperature goes to zero. 
Another quantity which is interesting also from the experimental point of view represents 
the specific heat. The thermal variations of this quantity are depicted in Fig. 6b for 
the same values of H/J as for entropy in Fig.6a. The displayed behaviour indicates a 
round Schottky-type maximum, whereas the stronger the external field, the flatter and 
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the broader the maximum. Apart from this trivial finding, one can also observe the 
double-peak specific heat curves (see Figs. 6c and 6d). The first peak which occurs in the 
specific heat curve at lower temperature is evidently closely related to the rapid variation 
of the magnetization (compare Fig.6c-d with Fig.5b). Moreover, it turns out that the 
maximum of this peak can be located approximately in the middle of ferrimagnetic region 
(in our case around H/ J m 1.5). On the other hand, the second peak may be thought as 
a Schottky-like peak resulting from the antiferromagnetic short-range order. Indeed, the 
relevant thermal dependencies for the nearest-neighbouring correlations strongly support 
this statement. 

In the following subsection we examine the spin S = 1 chain in order to clarify the 
influence of decorating spin on the magnetic properties of the system. 

3.2 Spin S = 1 chain 

We start our discussion once again with the analysis of ground state. In order to establish 
correct ground-state phase boundaries all relevant quantities have been examined in detail. 
From this analysis one can conclude that depending on the ratio between J 1; J (A) and 
H, altogether six different phases can appear in the ground state (see ground-state phase 
diagrams in Figs. 7a and 9a-b). In this part, we will firstly describe details of the spin 
ordering emerging in the appropriate ground-state phases and then, we will show how the 
scenario of the magnetization process depends on the parameters of the model. 

In the antiferromagnetic phase (AF) one finds a perfect antiferromagnetic alignment 
in the Ising sublattice (i.e. m s A = 1/2 and qff = —1/4) regardless of the strength of 
the exchange anisotropy A. Accordingly, the relevant spin order in the Ising sublat- 
tice is completely identical to that one in the AF phase of the spin S — 1/2 chain. 
Nevertheless, in contrast to the spin S — 1/2 case, the antiparallel alignment between 
nearest-neighbouring Heisenberg spins is weakened along z-axis as A increases from zero. 
In fact, the correlation qf^ tends monotonically from its classical Ising value q^ h = — 1 
at A = to q^l = —1/2 for large A (see e. g. Fig. 8a). In addition, one easily proves 
also the validity of relation rj = \qf^\ in the whole AF region. This observation would 
suggest that A supports the spin reorientation of Heisenberg spin pairs, namely, from 
the antiparallel oriented Heisenberg spin pair (one spin in S z = — 1 state, another one in 
spin S z = 1 state to be further denoted as the '1 — 1' spin pair) towards the spin pair 
with both Heisenberg spins in the spin S z = state (the '00' spin pair). However, in the 
large A limit both types of the Heisenberg spin pairs ('1 — V as well '00') are equally well 
populated and with a high probability also randomly distributed among Ising spins. This 
suggestion is strongly supported by results for the correlation qffi and staggered magne- 
tization m s B which asymptotically tend to zero as A — > oo (Fig. 8a). Finally, one should 
also notice that all these effects could originate from the onset of the antiferromagnetic 
short-range-ordering in the xy plane (nonzero qf^) that is the stronger (up to the value 
l/y/2), the greater the exchange anisotropy strength A. 

Now, let us turn our attention to both the FI and FII ferrimagnetic phases. As in the 
spin S — 1/2 case, the nonzero total magnetization of the FI phase arises due to the fully 
polarized Ising spins only {m z A = 1/2 and q zz = 1/4), while the Heisenberg sublattice 
does not contribute to the total magnetization at all (m z B = and qffi = 0). Since the FI 
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phase can arise by increasing the external field from the AF phase only (see Fig. 7a, 9ab), 
it is of interest to compare the relevant spin ordering in both phases. Actually, one still 
finds i] = \q^h\ to be valid, however, the nearest-neighbour correlation qfSi{<lhh) is weaker 
(stronger) in the FI phase with respect to that one in the AF phase (Fig.8d). These results 
are taken to mean that the number of nearest-neighbouring '00' spin pairs increases (of 
course, only up to one half of the total number of pairs) as one passes through the AF-FI 
phase boundary. Moreover, the appearance of a massive short-range order in the xy plane 
strongly implies the relevance of local quantum fluctuations also in the FI phase. On the 
other hand, in the second ferrimagnetic phase FII we have found the following results 

(C Qhh, q%, Qu, V, m z A , m' B , m) = (0, 1, -1/2, 1/4, 1, -1/2, 1, 1/2), 

indicating the 'classical' character of this phase that is usually observed also in the pure 
Ising spin systems. As the relevant spin-ordering is thoroughly analogous to that one in 
the FII phase of spin S = 1/2 chain, for brevity we will not repeat its description here. 

The most interesting spin order, however, can be found in the valence-bond phase 
(VB) and the intermediate phase (IP). Really, for instance in the VB phase one finds 

(qr h ,q z h z h ,q Z iH,Qu,V,rn z A ,m z B ,m) = (-1/2,0,1/4,1/4,1/2,1/2,1/2,1/2). 

As one can see, in contrast to the fully polarized Ising sublattice {m z A = 1/2 and qff = 
1/4), each couple of the nearest-neighbouring Heisenberg spins consists of one polarized 
spin (S z = 1) and one spin in the S z = state, i. e. the spin '01' pair. It is interesting 
to note that the symmetrization of both Heisenberg spin states can be achieved using 
the valence-bond-solid (VBS) picture [12]. Accordingly, each spin-1 atom splits into one 
polarized spin-1/2 variable with the fixed projection into the external field direction and 
one spin-1/2 variable with the unfixed projection creating a valence-bond. Thus, it is 
reasonable to assume that both Heisenberg spins interchange their spin states (tunneling 
between the S z = and S z = 1 spin states) and therefore, they effectively act on the 
surrounding Ising spins as the spin-1/2 atoms (see also the result for correlation qf%). 
However, the antiferromagnetic short-range order in the xy plane (nonzero <j^jQ, as well 
as the fact that the VB phase cannot be detected in the pure Ising system (i.e. for A = 0), 
imply again the obvious influence of the local quantum fluctuations. 

Perhaps the most striking spin alignment has been discovered in the IP phase. Among 
other matters, the perfect antiferromagnetic alignment between the Ising spins {m s A = 1/2 
and q zz = —1/4) has been confirmed in IP phase, hence, the total magnetization of system 
(m = 1/3) is entirely determined by the contribution of Heisenberg spins. Anyway, the 
results for the Heisenberg sublattice magnetization m z B = 1/2, correlation function qf^ = 
and quadrupolar momentum r\ = 1/2, clearly indicate that each couple of the nearest- 
neighbouring Heisenberg atoms comprises of the spin '01' pair. However, since both 
Heisenberg spins are placed in the IP phase between two non-equivalent Ising spins (one 
'up', other 'down'), the spin state interchange between both Heisenberg atoms would lead 
to qf% = what is in contradiction with our numerical result qffi ^ 0. On the other 
hand, if the Heisenberg atom in spin S z = 1 state would be strictly antiferromagnetically 
coupled to its nearest-neighbouring Ising spin (spin 'down'), then we would have = 
— 1/4. Nevertheless, our result for the correlation \q?fi\ < 1/4 indicates only the partial 
antiferromagnetic order between Ising and Heisenberg spins. This observation would 
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suggest that the spin '01' Heisenberg pairs should be, to a certain degree, randomly 
distributed among Ising spins. Finally, as one can expect, in the high field limit the system 
undergoes a phase transition towards the fully saturated paramagnetic phase (SP). As 
before, in the SP phase all Ising and Heisenberg spins are completely aligned towards the 
external-field direction, thus, in the SP phase one attains 

Qhl «£, C V, m' A , m%, m) = (0, 1, 1/2, 1/4, 1, 1/2, 1, 5/6); 

Now, let us proceed to examine the magnetization process of the system under inves- 
tigation. For this purpose, we have shown in Figs. 7a and 9a-b the ground-state phase 
diagrams in the A — H/ J plane together with some typical examples of the magnetization 
curves (Figs.7b-d and 9c-d) for three selected values of interaction parameters J\/J. In 
order to provide an independent check of the magnetization scenario from Figs.7b-d, the 
corresponding field-dependencies of the correlation functions and quadrupolar momentum 
are displayed in Figs.8b-d. Obviously, all the above results are absolutely in accordance 
with the aforementioned ground-state spin ordering. Moreover, through the comparison 
of Figs. 7a and 9a-b one can realize that the nature of magnetization process depends 
basically on the ratio between the Ising and Heisenberg interaction parameters. In fact, 
the stronger the Ising interaction J\ with respect to the Heisenberg one J (A), the broader 
the parameter region corresponding to the 'classical' FII phase. Otherwise, the increasing 
influence of the Heisenberg interaction J (A) causes the broadening of regions correspond- 
ing to the FI and VB phase, respectively, until the FII phase completely vanishes below 
Ji/J < 2/3, see e. g. Fig. 9a where the FII phase is already missing. Consequently, in the 
FI and VB phases one can expect that the effect of local quantum fluctuations plays an 
important role. 

Although, the system exhibits a stepwise magnetization curve with an abrupt change 
of the magnetization in the whole range of parameters (see Figs.7b-d and 9c-d), there is 
a fundamental difference between the magnetization process in the 'classical' Ising-like 
regime and that one in the quantual regime. As a matter of fact, in the former case one 
observes the double-plateau magnetization curve with the FII plateau state only, i. e. 
with the FII phase as an intermediate state between the AF and SP phase (see Fig. 7b). 
Contrary to this, in the latter case one can detect the double- (Fig. 7c), triple- (Fig.7d), 
or quadruple-plateau (Fig.9c-d) magnetization curves. The double-plateau magnetiza- 
tion curves are rather rarely observable, since they arise due to the direct metamagnetic 
transition from the AF phase to the VB phase in a relatively narrow region of A only 
(see Fig. 7a). However, when the mechanism of the magnetization process is driven by 
the Heisenberg interaction, i. e. under the requirement of sufficiently small J\j J and 
sufficiently large A, the triple-plateau magnetization curves with the transitions between 
AF-FI-VB-SP phases are always preferred (Fig.7d). Finally, one should also remark that 
if the condition A m H/J m Ji/J> 1 is satisfied, there also appear the extraordinary 
quadruple-plateau magnetization curves (see Fig. 9b) with the IP phase in a very narrow 
region of the external fields. Indeed, we found even two different possibilities for the 
quadruple-cascade transitions, namely, the AF-FI-IP-FII-SP cascade transitions (Fig. 9c) 
as well as the AF-FI-IP-VB-SP ones (Fig.9d). 

To conclude the analysis of the spin S = 1 chain, we will also briefly mention the 
finite-temperature behaviour of the system under consideration. For this purpose, the 
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thermal variations of the total magnetization are plotted in Fig. 10 for Ji/J = 1.0 and two 
selected values of the anisotropy A. In agreement with the aforementioned arguments, 
one observes here three and two field-induced transitions for the anisotropy strengths 
A = 1.0 and 0.55, respectively. Evidently, as the magnitude of external field varies, the 
various thermal dependencies result from the competition between the Ising interaction, 
Heisenberg interaction and magnetic field. However, the most interesting dependencies 
arise again for the external fields from the vicinity of the phase boundaries. In such a 
case, the relevant thermal excitations result in a very rapid change of magnetization due 
to the competing influence of the phases separated by relevant transition line. Moreover, 
it turns out that the narrower interval of external fields corresponds to the relevant phase, 
the more robust change in the magnetization can be observed. 

4 Conclusion 

In the present article we have obtained the exact solution of the mixed spin-1/2 and spin- 
s' (S = 1/2, 1) Ising-Heisenberg chain in the external magnetic field. The most important 
result stemming from this study is the confirmation of multistep magnetization process 
by an exact calculation. Moreover, it has been proved that the character of magnetization 
process depends essentially on the ratio between the Ising and Heisenberg interactions, 
whereas the XXZ anisotropy term A also allows to control it in deciding manner. Since 
the presence of the non-diagonal interaction term J A is responsible for the onset of the 
local quantum fluctuations, as we have also shown, it basically modifies an otherwise 
trivial Ising-like behaviour. For example, in the case of the spin S = 1 chain one finds 
instead of the double-plateau magnetization curve arising in the pure Ising spin system 
(A = 0), the double-, triple-, or even quadruple-plateau magnetization curves in the 
Ising-Heisenberg chain with A ^ 0. Altogether, the presented results indicate that the 
extraordinary rich ground-state phase diagrams result from the competition between the 
easy-axis interactions J±, J and the easy-plane interaction J A. 

One should also emphasize that our research on Ising-Heisenberg chains has been stim- 
ulated by the recent experimental works dealing with many quasi-lD mixed-spin chains 
[13]. Although, we are not aware of any quasi-lD system in that two kinds of magnetic 
ions regularly alternate with p = 3 periodic fashion (AABAAB...), the recent progress in 
molecular engineering [14] supports our hope that the synthesis of such a polymeric chain 
should be possible in the near future. Structural derivatives of a novel polymeric chain 
recently reported by Mukherjeee et al [15], seem to be the most promising candidates 
from this point of view. In fact, the crystal structure of the above mentioned polymeric 
chain consists of the spin-1/2 Cu 11 dimers linked through Ni 11 monomers. Unfortunately, 
as a consequence of the square-planar coordination of Ni ions in [Ni n (CN) 4 ] 2 - bridg- 
ing groups, the Ni 11 metal ions are diamagnetic and thus, they do not contribute to the 
magnetism. 

Finally, it should be stressed that the applied mapping transformation technique does 
not require any restriction to the dimensionality of the spin system and hence, it can be 
straightforwardly generalized to the mixed Ising-Heisenberg lattices in two- and three- 
dimensions, too [10]. 
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5 Appendix 

Explicit expressions for the functions V\-V G . 

a) Spin- 1/2 chain: 

Vx = 2 exp(-/3 J/4) cosh(/3 Ji/2 + (3H B ) + 2 exp(/3J/4) cosh(/3 JA/2) 
V 2 = 2exp(-/3J/4) cosh(/3Ji/2 - (3H B ) + 2exp(/3J/4) cosh(/3JA/2) 

V 3 = 2exp(-/3J/4) cosh(/3/J B ) + 2exp(/3J/4) cosh(/^J 2 + (JA) 2 /2) 
V 4 = exp(-/3 J/4) sinh(/5Ji/2 + /3/f B )/2 
\/ 5 = exp(-/3J/4) sinh(/3 Ji/2 - (3H B )/2 
V 6 = exp(-/3J/4) sinh(/3JT B )/2, 

b) Spin-1 chain: 

Vi = 2 exp(-/3 J) cosh(/3 Jx + 2/3/J B ) + exp(2/3 J/3) Wi + 

+4 cosh(/3 Jx/2 + /3# B ) cosh(/3 JA) 
F 2 = 2 exp(-/3 J) cosh(/3 J x - 2(3H B ) + exp(2(3J/3)W 1 + 

+4cosh(/3Ji/2 - (3H B ) cosh^JA) 
\/ 3 = 2exp(-/3J) cosh(2/3/J B ) + exp(2/3J/3)W 2 + 

+4cosh(/3/J B ) cosh(/3y/j? + (2JA) 2 /2) 
Vk = exp(-/3J) sinh(/3 J x + 2/3if B ) + sinh((3J 1 /2 + /3if B ) cosh^JA) 
V 5 = exp(-/3J) sinh(/3Ji - 2/3if B ) + sinh(/?Ji/2 - /3if B ) cosh(/3JA) 
\/ 6 = exp(-/3J) sinh(2/3/J B ) + sinh(/3/J B ) cosh^Jf + (2JA) 2 /2) 

where the expressions W\ and W 2 are given by: 

2 

Wi = exp{-2/?Pi cos[(0i + 27rn)/3]}, 

n=0 
2 

W 2 = J2 exp{-2/3P 2 cos[(0 2 + 2vrn)/3]}, 

n=0 

Pi = (J/3) 2 + 2( JA) 2 /3, P| = (J/3) 2 + 2( JA) 2 /3 + 

Qi = (J/3) 3 + J(JA) 2 /3, Q 2 = (J/3) 3 + J(JA) 2 /3 - J 3 , 

0x = arctan (yJP?-Ql/Qi) , 2 = arctan ( ^ P$ - Q 2 2 /Q 2 ) ■ 
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Figure captions 

Fig.l Part of the doubly decorated mixed-spin chain. The black circles denote the spin- 
1/2 Ising atoms of sublattice A and the gray ones represent the decorating spin-S' 
Heisenberg atoms of sublattice B. The ellipse demarcates a typical bond described 
by the Hamiltonian Hk introduced in Eq. (2). 

Fig. 2 Ground-state phase boundaries in the A - H/ J plane for Ji/J = 1.0 and 2.0. 

Fig. 3 Low-temperature (ksT/J = 0.01) magnetization curves for: a) Ji/J = 1.0 and 
A = 1.0; b)-d) J^J = 2.0 and A = 0.5,1.0,1.5. The solid and dotted (dashed) 
lines represent the total magnetization per one site and the single-site magnetization 
of the Ising (Heisenberg) sublattice, respectively. 

Fig. 4 Low-temperature behaviour [k B T j J = 0.01) of several correlation functions: qff h 
and qff (dashed lines), qf^ (dotted lines) and qffc (solid lines). 

a) plot of the correlation functions against the anisotropy A for the system without 
external field; 

b) -d) field-dependence of the correlation functions depicted for Ji/J = 2.0 and 
various anisotropies A = 0.5, 1.0 and 1.5. 

Fig. 5 Thermal variations of the total magnetization for some typical values of the external 
field, Ji/J = 2.0 and two selected anisotropies A = 0.5 and 1.5, respectively. 

Fig. 6 Thermal behaviour of the system for Ji/J = 2.0, A = 1.5 and several values of 
external field. 

a) entropy plot versus temperature; 

b) -d) variations of the specific heat with temperature. 

Fig. 7 Ground-state phase diagram together with some typical examples of low-temperature 
magnetization curves when the ratio J\/J— 1.0; 

a) ground-state phase diagram in the A — H/ J plane; 

b) -d) some typical examples of the low-temperature (/c^T/J = 0.001) magnetization 
curves for various anisotropies A. 

Fig. 8 Low-temperature (ksT/J = 0.001) behaviour of several correlation functions qf^ 
(dashed), qf^ (dotted), q^ (solid), qff (dashed-dotted), quadrupolar momentum rj 
(dashed-dotted-dotted) and staggered sublattice magnetization m s B (solid line). 

a) zero-field variations of relevant quantities with anisotropy A; 

b) -d) field-dependencies of correlations for selected values of A. 

Fig. 9 Ground-state phase diagrams together with some typical examples of low-temperature 
magnetization curves: 

a)-b) ground-state phase diagram in the A — H/ J plane for J\j J = 0.5 and 2.0; 

c) -d) some typical examples of the low-temperature (ksT / J = 0.001) magnetization 
curves when Ji/J = 2.0 and A = 1.9 and 2.05, respectively. 

Fig. 10 Thermal variations of the total magnetization for J\j J = 1.0 and two selected values 
of A. 
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